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Abstract 

We study the dihadron azimuthal correlation produced nearly back-to-back in unpolarized 
hadron collisions, arising from the product of two Collins fragmentation functions. Using the 
latest Collins fragmentation functions extracted from the global analysis of available experimental 
data, we make predictions for the azimuthal correlation of two-pion production in pp collisions at 
RHIC energies. We find that the correlation is sizable in the mid-rapidity region for moderate jet 
transverse momentum. 
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I. INTRODUCTION 



The transverse momentum dependent (TMD) distributions and fragmentation functions 
have received much attention recently jl|. They are believed to be responsible for several 
azimuthal asymmetries observed in experiments, such as the single transverse spin asym 
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metry (SSA) in semi-inclusive hadron production in deep inelastic scattering (SIDIS) [2|, 
and in hadronic collisions [4], as well as the large cos(20) anomalous azimuthal asymmetry 
in back-to-back dihadron production in e + e~ annihilation J5]. 

Among these TMD parton distributions and fragmentation functions, the Sivers quark 

fl Q 

distribution |6| and the Collins fragmentation function |7| are mostly discussed in the last 

few years. The Sivers quark distribution represents a distribution of unpolarized quarks 
in a transversely polarized nucleon, through a correlation between the quark's transverse 
momentum and the nucleon polarization vector. On the other hand, the Collins fragmenta- 
tion function describes a transversely polarized quark jet fragmenting into an unpolarized 
hadron, whose transverse momentum relative to the jet axis correlates with the transverse 
polarization vector of the fragmenting quark. 

Though both of them belong to the so-called "naive-time-reversal-odd" (T-odd) func- 
tions, they have very different universality properties. For the Sivers functions, it has been 
shown that they differ by a sign for the SIDIS and Drell-Yan (DY) processes [8] , and those 
in the hadronic collisions have even more nontrivial relation to that in SIDIS and DY pro- 
cesses 9Nl3|. On the other hand, the Collins fragmentation function is universal between 
different processes, in the SIDIS, e + e~ and hadronic collisions 14l-ll6l]. The effect of the 
Collins fragmentation function has been recently explored by one of us in the azimuthal 
asymmetric distribution of hadrons inside a jet in p^p collision [l^]. It is demonstrated that 
the asymmetry is sizable at RHIC, therefore, the experimental study of this process could 
provide an important information on the universality of the Collins fragmentation function. 

Another difference between these two functions is that the Collins fragmentation function 
is chiral-odd whereas the Sivers function is chiral-even. Because of its chiral-odd nature, the 



Collins effect can only be observed when it is coupled to another chira i 
or fragmentation function. In SIDIS, the chiral-odd quark transversity 
the Collins fragmentation function and leads to nonzero azimuthal SSA 



-odd distribution 
can couple to 
7fl. This SSA has 

been studied by the HERMES || and COMPASS |3| collaborations, and very interesting 
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results on the Collins fragmentation function have been found. In e + e~ annihilation pro- 
cess, two Collins fragmentation functions couple to each other in the back-to-back dihadron 
production, results into to a cos(2</>) azimuthal asymmetry jl9 |. This anomalous cos(2</>) 
asymmetry has been measured by the BELLE Collaboration 5[, and was found consistent 
with the HERMES and COMPASS measurements on the Collins fragmentation functions. 



Recently a global analysis of these experimenta 



fragmentation functions have been extracted 



20] 



data has been performed and the Collins 



In this paper we investigate the possibility of exploring the Collins fragmentation function 
in unpolarized pp collision by studying the azimuthal correlation in back-to-back dihadron 
production, following the same wisdom of dihadron production in e + e~ annihilation. We 
show that the asymmetry is proportional to the product of two Collins fragmentation func- 
tion, same as that in e + e~ annihilation. Using the latest Collins fragmentation function 
extracted from the global analysis of available data on SIDIS and e + e~ experiments, we es- 
timate the asymmetry for dihadron production at RHIC energy. We find that the azimuthal 
asymmetry is sizable at mid-rapidity region for moderate jet transverse momentum. We 
argue that this process shall provide addtional important information on the Collins frag- 
mentation function and its universality properties. 

The rest of the paper is organized as follows. In Sec. [TTJ, we derive the theoretical results 
for the dihadron azimuthal correlation produced nearly back-to-back in unpolarized hadron 
collision. In Sec. IHIt we present our numerical predictions for the azimuthal correlation in 
unpolarized pp collisions for RHIC kinematics. Finally, we summarize our findings and the 
corresponding conclusions in Sec. [TV] 



II. DIHADRON AZIMUTHAL CORRELATION IN UNPOLARIZED HADRON 
COLLISION 

We study the azimuthal correlation of two hadrons h± and h 2 produced nearly back-to- 
back in a hadronic collision, 

A(P 1 ) + B(P 2 ) h(P hl ) + h 2 (P h2 ) + X, (1) 

where both of the incident hadrons A and B are unpolarized. The momenta of the initial 
hadrons are denoted by Pi and P 2 , and those of the final hadrons by Phi and Ph 2 , respectively. 
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FIG. 1: The leading order contribution to the cross section of A(P\) + B(P2) — > hi(Phi) + ^2(^2) + 
X with the 2 — > 2 partonic process a(k±) + 6(^2) — > c{q\) + <i(q 2 ). 

The leading order contribution to the scattering cross section comes from partonic 2—7-2 

sub-processes, a{ki) + b(k 2 ) — > c(q\) + <i(g 2 ), as shown in Fig. [TJ The parton momenta are 
expanded as follows, 

k x = x x P\ + k 1T , (2a) 

A; 2 = x 2 P 2 + fc 2T , (2b) 

Phi = ziqi + Pit, (2c) 

Ph2 = z 2 q2 + P2T, (2d) 

where X\ and x 2 are the longitudinal momentum fractions, and ky? and /c 2 t are the transverse 
momentum of the parton relative to the corresponding incident hadron. qi and g 2 are the 
momenta of the nearly back-to-back jets J\ and J 2 , which has a polar angle 9\ and 6* 2 relative 
to the incoming hadron Pi, respectively. The momenta of the incoming hadrons and the 
final state two jets form the so-called reaction plane (approximately). Besides carrying a 
longitudinal momentum fraction Z\ (z 2 ) of the jet Ji (J 2 ), the hadron hi (h 2 ) also has 
a transverse momentum pix (p 2 t) relative to the jet J\ (J 2 ) direction, which defines an 
azimuthal angle with the reaction plane: <pi (0 2 ), as shown in Fig. |2J Due to the Collins 
effect, there will be a correlation between these two azimuthal angles <p\ and 2 , which 
is proportional to the product of the two Collins fragmentation functions as we will show 
below. 
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FIG. 2: Illustration of the kinematics for dihadron production AB — > h\h 2 + X. 



From Fig. [T], the cross section of dihadron production can be written as [l3[ , 



da 



dx\d 2 kiTdx 2 d 2 k^dzx d 2 pixdz 2 d 2 p 2 r - 



d3qi d3q2 -WS'ih 



2S j _ rx™- (27r)32jEl (27r) 3 2E 2 v "' v " ^ ' ^ 91 ^ 

xTr /ci T )$(x 2 , k 2T )A(zi,pi T )A(z 2 ,p 2T )H(ki, k 2 , q u q 2 )H*(k 1 , k 2 , q u q 2 )] , (3) 



where S = (Pi + P 2 ) 2 , and &t) and A(z,pt) are the distribution and fragmentation 
correlation functions, and H(ki,k 2 ,qi,q 2 ) are the hard part amplitude. 

? or the unpolarized hadron, the correlation functions hp) can be simply decomposed 



as 
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22] 



$(P,z,fcr) 



f( Z ,k 2 T )r+K{z,k 2 T ) 



2 .a^k T ^P u 



M 



(4) 



where f(x,kj,) is the unpolarized TMD parton distribution function, and hi( the 



Boer- Mulders function |22j . Similarly, we can parameterize the gluon distributions from the 
incoming hadrons. The effect of Boer- Mulders function in unpolarized hadronic collisions has 
been extensively studied previously [23|-|25| , which will be neglected in our current study. We 
will concentrate on the effect coming from the fragmentation correlation function A(z,pt), 
which can be expanded as [26] 



A(P h ,z,p T ) 



D(z,p 2 T )P h + Ht(z,p 2 T ) 



2 ^ v P hil p Tv 



zM h 



(5) 



where D(z,pj,) is the unpolarized TMD fragmentation function, and is the Collins 
function. From its definition, we can see that the Collins function describes a transversely 
polarized quark jet fragmenting into an unpolarized hadron jj]. It is this Collins function 
that generates a non-vanishing azimuthal correlation between the final state two hadrons. 
Since Collins function is a chiral-odd TMD function, the azimuthal correlation will depend 
on the product of two Collins functions. 
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The phase space integral in Eq. fl3]) can be simplified by using 
d 3 qi d 3 q 2 



(2n) 3 2E 1 (2tt) 3 2£ 2 
1 



dy 1 dy 2 dP 1 



(2vr) 4 -5\k 1 + k 2 -q 1 -q 2 ) 
P 



8nS 



s 



Vs 



-ML 



+ e 



in - 



s 



x 2 



Pj - ( - 



yi 



(6) 



where y\ and y 2 are rapidities for the jet J\ and J 2 , P± is the jet transverse momentum. 
Since Pit,P2T *C P±, the rapidity of the hadron is approximately equal to that of the parent 
jet. As stated above, we neglect the intrinsic transverse momentum effects for the incoming 
partons. Finally we obtain the cross section for this process as 

da 



dVS 



Efa{%l) fb{x 2 ) r 
— — [Dc^hAZliPlT) D d^hA z 2,P 2 T) H ab^cd \S,t,u) 



abed 

I Pit ■ P2T 



x 2 

X\X 2 

~pj 



(Pi • PytPi ■ P2T + Pi ■ P2TP2 ■ Pit) 



Hi(Zi,pI t )Hi(z 2 ,pI t ) ^Collins (t, 1 
X — — h ab^cd \S,t,U)\ 



(7) 



Zl M h z 2 M h 

where dVS = dyidy 2 dP]_dzidz 2 d 2 pird 2 p 2 T is the phase space for this process, f a (xi) and 
fb{x 2 ) are the standard unpolarized parton distribution functions, and s, t, and u are the 
usual partonic Mandelstam variables. The parton momentum fractions X\ and x 2 are fixed 
by the delta functions in Eq. (jHJ), 



x 2 



P± 

Vs 
p± 
Vs 



(8a) 
(8b) 
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The normal partonic cross sections H^ b ^ cd are well-known 27J] 
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(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



The new hard parts H%gj%g that are responsible for the azimuthal correlation are given by, 



H 



Collins 

gg^qq 



N c N 2 - 1 



2tu 



H, 



Collins 

qq-^q'q' 



N 2 
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32 



N 2 



N 3 

c 



H 



Collins 

qq->qq 



0, 
N 2 



(17) 

(18) 

(19) 
(20) 
(21) 



and the hard parts for the channels with gluon in the final state, qq — > gg, qg — > qg and 
99 99 vanish since there is no gluon Collins function. 

III. PHENOMENOLOGY STUDY 



In this section we first properly define the azimuthal asymmetry to be measured in the 
experiments. We then use the latest Collins fragmentation function to estimate this asym- 
metry for RHIC kinematics. 
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From Eq. <^j, the explicit form of the azimuthal correlation depends on the following 
function, 

h(pi T ,P2T, 0i, 02) = Pit • Pit - • P1T-P2 • P2T + Pi ■ P2TP2 • Pit)- (22) 

M. 

Function h(piT,P2T, <f>i, $2) is reduced to its simplest form in the partonic center of mass 
frame, 

KP1T,P2T, 01, 02) = -P1TP2T cos(0i + 2 ) . (23) 

To take advantage of this simplicity, one could boost from the Lab frame to the partonic 
CM frame experimentally. Or equivalently, one can select the events with y\ + y 2 ~ 0, where 
the Lab frame coincides with the partonic CM frame. In the rest of our paper, we will 
take yi + y 2 = in our calculations and present the numerical estimate for the azimuthal 
asymmetry. 

Integrating over the intrinsic transverse momentum p^T an d P2T1 we have 

do- vra 2 ^ f a (xx) f b {x 2 ) . , , , v * ~\ foA s 

d yi dy 2 dPld Zl dz 2 = ^ X * Dc ^ iZl)Dd ~* MH ^ cd (S ' M) (24) 

While integrating over the moduli of the intrinsic momenta pyt and p 2 T, and over the 
azimuthal angle 0i, one obtain 

dy l dy 2 dPldz l dz 2 d^ l + ^ 2 ) = 2S^^~Ti —L^Ui)^(^^ 

-oos^ + ^^W^C^SS^*.")] . ( 25 ) 
where ^g'- 1 / 2 -'^) is the so-called half-moment of the Collins function given by, 

} (i/2) m _ /" j2_ „ H i(z,Pt) 



6q^(z) = I d 2 PTPT vy zM ^' . (26) 



Then following the normal analysis in e + e — >■ /ii/i 2 + X, we define 

(yi,y 2 ,P±, Z U Z 2 , 0! + 2 ) 

da / 1 da 



dy\dy 2 dP\dz\dz 2 d (0i + 2 ) / 27r dy\dy 2 dP\dz\dz 2 



s 



Again to eliminate the false asymmetries we introduce the ratio of the unlike-sign to 
like-sign pion production, Ajj and Al, given by 

1 - cos(0i + 4> 2 )Pu 



R 



Au 
A L 
1 - 



1 - COS(0! + 4>2)Pl 
COS(01 + 02)^12(2/1, 2/2, P±, Zl, Z 2 ) 



1 - COS(0 1 + (j> 2 ) (P V - P L ) 



(28) 



with 



v 







rrCollins 
abated 


J2abcdfa(Xl)fb(x 2 ) 
T,abcdfa(Xl)fb(x2) 


[D c ^+(z 1 )D d ^ 7r -(z 2 ) + D c ^ n -(z 1 )D d ^ w +(z 2 )} 


abated 

rrCollins 
abated 


J2abcdfa{xi)fb(x 2 ) [D c ^+{zi)D d ^+(z 2 ) + D c ^-{z 1 )D d ^-{z 2 )\ H^ cd 



Pi 



-4l2(2/l, 2/2, PjL, Zl, Z 2 ) = - ^L- 



(29) 

(30) 
(31) 



This way the true asymmetry due to Collins effect is encoded in the so-called double 
ratio asymmetry parameter Ai 2 (yx, y 2 , P±, z±, z 2 ). To evaluate A 12 for the dihadron produc- 
tion in un pol arized pp collision at y/S = 200GeV at RHIC, we use the Collins fragmentation 

ERMES, COM- 
2811 and Kretzer 



functions [20[ extracted from a combined fit to the experimental data from £ 
PASS and BELLE collaborations. We use CTEQ5L parton distributions 



unpolarized fragmentation function obtained in 

cm 0.05 



29(]. 



0.04 
0.03 
0.02 
0.01 



P ± =4 GeV 



, i , , , , i , , . , i . , , , i , 




, I , , , , I , , , . I , 



0.05 
0.04 
0.03 
0.02 
0.01 



P ± =6 GeV 



, i , , , , i , , . , i . , , , i , 




, I , , , , I , , , . I , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 3: Azimuthal asymmetry ratios Ai 2 (yi, y 2 , P±, z±, z 2 ) (defined in Eq. (|28|) ) of unlike-sign to 
like-sign pion production in unpolarized proton-proton scattering at RHIC energy y/~S = 200GeV, 
as functions of z 2 with 2/1 = 2/2 = 0, P± = 4 GeV (left) or P± = 6 GeV (right). The curves are: 
solid (0.7 < zi < 0.9), dashed (0.5 < z 1 < 0.7) and dotted (0.3 < z 1 < 0.5). 
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In Fig. |3l we plot Ai 2 (yi,y 2 , P±, Z\, z 2 ) for the dihadron production in mid-rapidity y\ = 
y 2 = region at RHIC energy y/S = 200 GeV. In the figure on the left, we plot A i2 as a 
function of z 2 at the jet transverse momentum Pj_ = 4 GeV with Z\ integrated from three 
different ranges 0.3 < z\ < 0.5 (dotted), 0.5 < z\ < 0.7 (dashed), and 0.7 < z\ < 0.9 (solid), 
respectively. On the right, we present the same plot but with P± — 6 GeV. We find that the 
asymmetry A± 2 is largest when both Z\ and z 2 become large, same as what has been observed 
in e + e~ experiments [5J. On the other hand, A 12 decreases when increasing P±. This is also 
consistent with what BELLE observed if one realizes that s/APj_ = sin 2 9 in parton CM 
frame. Though it has similar features as that in e + e~ collision, the asymmetry in hadronic 
collision is smaller. This is due to the fact that there is copious gg — > gg and qg — > qg 
contribution to the azimuthal angle independent cross section, while they do not contribute 
to the azimuthal dependent part since there is no gluon Collins function. However, the 
asymmetry is still around several percent and shall be measurable at RHIC. 

These results can be extended to general kinematics, for example, in two different rapidity 
regions: \yi\ ^ \y 2 \. in this case, although the azimuthal angular dependence is not exactly 
as cos(0i + 4> 2 ) in Eq. ( l23i) . the Collins fragmentation functions will nevertheless lead to a 
nonzero mean value of (cos(0i + <f) 2 ))- This can be seen from the differential cross section 
expression in Eq. ([7]). The experimental observation of this nonzero effects can be used 
as signal of the Collins effects, since the normal fragmentation functions D(z,px) will not 
contribute to a nonzero (cos(0i + 4> 2 ))- We hope that the future RHIC experiments can 
carry out these measurements, and provide more information on the Collins fragmentation 
functions, which will help us to pin down the mechanism for the single spin asymmetry in 
hadronic collisions as we discussed in the Introduction. 



IV. SUMMARY 



In this paper, we have studied the dihadron azimuthal correlation produced nearly back- 
to-back in unpolarized hadron collision, arising from the product of two Collins fragmen- 
tation functions. Using the latest Collins fragmentation function extracted from the global 
analysis of available experimental data, we make predictions for the azimuthal correlation 
of two-pion production in unpolarized pp collisions at RHIC energies. We find that the 
feature of the asymmetry is similar to those observed in e + e~ annihilation. The asymmetry 
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parameter A\2 is sizable in the mid-rapidity region for moderate jet transverse momentum, 
and could be measurable in the experiments. The experimental study of this process could 
provide the important information on the size of Collins fragmentation function in hadronic 
collision, at the same time, it could also be used to test the universality properties of the 
Collins fragmentation function in different processes. 
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